Ferric iron partitioning between plagioclase and silicate liquid: thermodynamics and petrological applications

A series of Fe and Mg partition experiments between plagioclase and silicate liquid were performed in the system SiO₂-Al₂O₃-Fe₂O₃-FeO-MgO-CaO-Na₂O under oxygen fugacities from below the IW buffer up to that of air. A thermodynamic model of plagioclase solid solution for the (CaAl,NaSi,KSi)(Fe³⁺,Al³⁺)Si₂O₈-Ca(Fe²⁺,Mg)Si₃O₈ system is proposed and is calibrated by regression analysis based on new and previously reported experimental data of Fe and Mg partitioning between plagioclase and silicate liquid, and reported thermodynamic properties of end members, ternary feldspar and silicate liquid. Using the derived thermodynamic model, FeO^t, MgO content and Mg/(Fe^t+Mg) in plagioclase can be predicted from liquid composition with standard deviations of ǂ.34 wt% (relative error =9%) and ǂ.08 wt% (14%) and ǂ.7 (8%) respectively. Calculated Fe³⁺-Al exchange chemical potentials of plagioclase, m_{Fe^{3 +} ( Al )- 1} ^Pl{\rm \mu }_{{\rm Fe}^{{\rm 3 + }} \left( {{\rm Al}} \right)_{{\rm - 1}} }^{{\rm Pl}} agree with those calculated using reported thermodynamic models for multicomponent spinel, m_{Fe^{3 +} ( Al )- 1} ^Sp{\rm \mu }_{{\rm Fe}^{{\rm 3 + }} \left( {{\rm Al}} \right)_{{\rm - 1}} }^{{\rm Sp}} and clinopyroxene, m_{Fe^{3 +} ( Al )- 1} ^Cpx{\rm \mu }_{{\rm Fe}^{{\rm 3 + }} \left( {{\rm Al}} \right)_{{\rm - 1}} }^{{\rm Cpx}} . The FeO^t content of plagioclase coexisting with spinel or clinopyroxene is affected by Fe³⁺/(Fe³⁺+Al) and Mg/(Fe+Mg) of spinel or clinopyroxene and temperature, while it is independent of the anorthite content of plagioclase. Three oxygen barometers based on the proposed model are investigated. Although the oxygen fugacities predicted by the plagioclase-liquid oxygen barometer are scattered, this study found that plagioclase-spinel-clinopyroxene-oxygen and plagioclase-olivine-oxygen equilibria can be used as practical oxygen barometers. As a petrological application, prediction of plagioclase composition and f_O2 are carried out for the Upper Zone of the Skaergaard intrusion. The estimated oxygen fugacities are well below QFM buffer and consistent with the estimation of oxidization states in previous studies.     

Dislocation microstructures in naturally deformed silicate garnets

Dislocation microstructures of naturally deformed silicate garnets and olivines in garnet-peridotites and silicate garnets in eclogites from four localities have been observed with a transmission electron microscope (TEM) to clarify the dislocation characteristics of silicate garnets. We have obtained the following results: (1) dislocation densities of garnets in all the garnet-peridotites (ρ = 10⁵−10⁷ cm⁻²) are always nearly an order of magnitude lower than those of co-existing olivines; (2) dislocation densities of garnets in eclogites (ρ = 10⁵−10⁸ cm⁻²) which are embedded in garnet-peridotites are almost an order of magnitude higher than those of garnets in the surrounding garnet-peridotites; (3) the dominant Burgers vector, b, of mainly edge dislocations in garnet is 100 for specimens with dislocation density ρ = 10⁵−10⁶ cm⁻², while b=1/2111 for specimens with ρ = 10⁷−10⁸ cm⁻². Result (1) indicates that the observed dislocations in garnets were formed by plastic deformation under the same stresses as for co-existing olivines, and that there is a similar relationship between applied stress and dislocation density for garnets as for olivines. Result (2) suggests that the stress concentration occurred around eclogites embedded in garnet-peridotites, and the resulting differential stress in garnets in eclogites was further elevated by the interlocking of neighboring hard garnet grains. Finally, result (3) indicates that the dominant Burgers vector of mainly edge dislocations in garnet changes from 100 to 1/2111 with increasing applied differential stress.     

Relationship between faults and generation of inland shallow earthquakes: enlightenment of the M=7.2 southern Hyogo prefecture earthquake

The M=7.2 southern Hyogo prefecture earthquake, which occurred on January 17, 1995, destroyed the region severely. Further researches are necessary to explain the problems obtained from this earthquake. We have discussed some characteristics correlating with this earthquake, such as distribution of aftershocks and disaster, relationship between fault and generation of inland shallow earthquakes, strong motion and so on. The tectonic movements of surface faults may be different from those of deep faults. Earthquakes may not be restricted directly by the tectonic movements of surface faults. The strong motion zone is often appeared in the region a little away from the both ends or from the single side of a fault instead of the region of a focal fault. Some mechanisms of strong motion have also been discussed.     

Holocene beach deposits for assessing coastal uplift of the northeastern Boso Peninsula, Pacific coast of Japan

This paper presents a case study that assessed spatial variations in the tectonic uplift rates of beach deposits in the relict Kujukuri strand plain, situated on the northeastern coast of the Boso Peninsula, eastern Japan. The southern Boso Peninsula is tilted downward to the northeast due to plate subduction along the Sagami Trough. However, the cause of the northeastern coast uplift creating the relict strand plain is unclear, due to the absence of a Holocene raised marine terrace sequence. Elevations and ages of beach deposits were collected from drilled cores and ground-penetrating radar profiles along three shore-normal sections in the southern Kujukuri strand plain. From this, alongshore variations in the relative sea level since the mid-Holocene could be seen. These corresponded to north-to-northeast downward tilting at a rate of 0.4 m/ka for an interval 10 km and are concordant with the longer term tilting of the last interglacial marine terrace surrounding the plain. Although it is difficult to assess shore-normal variations of uplift based on the present dataset, the recognized tilting apparently continues to the tilting of the southern Boso Peninsula, implying the Sagami Trough probably affects the uplift of the Kujukuri coast.     

Atoll island vulnerability to flooding and inundation revealed by historical reconstruction: Fongafale Islet, Funafuti Atoll, Tuvalu

The reef islands formed on coral atolls are generally small, low, and flat, with elevations of only a few meters. These islands are thus highly vulnerable to elevated sea levels caused by extreme events and global warming. Such vulnerability was recently evidenced at Fongafale Islet, the capital of Tuvalu, when it flooded during accelerated spring high tides possibly related to sea level rise caused by global warming. Many factors, not only environmental but also economic and social, determine the vulnerability of an island to sea level rise. In this study, we used data spanning 108 yrs to reconstruct changes in topography, land use/cover, population, and the distribution of buildings at Fongafale Islet. The results indicate that the vulnerability of Fongafale Islet relates to its original landform characteristics: the central part of the island was formerly dominated by swampland that flooded at high tides. Fongafale Islet experienced greater population in-migration and centralization beginning in the 1970s following the independence of Tuvalu and Kiribati. Migrants were also responding to declines in overseas mining operations and limited options for paid employment. As the population increased, construction took place in vulnerable swampland areas. Our results clearly demonstrate that examinations of global environmental issues should focus on characteristics specific to the region of interest. These characteristics should be specified using historical reconstruction to understand and address the vulnerability of an area to global environmental changes.     

Collapse of primordial clouds and the formation of a central core

The collapse of rotating clouds is investigated using three-dimensional self-gravitating hydrodynamical simulations. We take into account the detailed non-equilibrium chemical reactions for primordial gas that consists of pure hydrogen. The parameters of the collapse and the condition of the fragmentation are compared with those of isothermal clouds. It is shown that the geometrical flatness of the central region of the disc is a good indicator for predicting whether the clouds fragment or not. If the flatness is greater than the critical value, ∼ 4π, a cloud fragments into filaments and blobs. On the other hand, if the flatness is smaller than the critical value, fragmentation is not expected before the central core formation even if the cooling is efficient and the total mass becomes much greater than the local Jeans mass at the centre. The critical mass is found to be 3 × 10⁶ M_⊙ for a typical initial condition. If the initial cloud mass is smaller than this critical value, fragmentation before the central core formation is not expected. For a typically estimated first collapsing cosmological baryonic object, M ≲ 10⁶ M_⊙, central core formation is expected and will have a significant effect on the later evolution of the whole system. This revised version was published online in July 2006 with corrections to the Cover Date.     

The eccentric-orbit binary model for the transient X-ray sources

Eccentric-orbit binary models for transient X-ray sources are investigated. In these models, a compact star is in an eccentric orbit around a more massive star. As the compact star accretes mass from the stellar wind of the massive star, the accretion rate becomes time-dependent. The accretion rate is determined by Bondi's accretion radius, which depends on both the relative velocity of the stellar wind to the compact star and the sound velocity through the stellar wind. With reasonable sets of the eccentricity, the semi-major axis, the stellar wind velocity and the sound velocity, we obtain the variations of the light curves compatible with observations for the transient X-ray sources. It is likely that many transient X-ray sources are explainable by eccentric-orbit binary models.     

Radiative Heat Transfer and Hydrostatic Stability in Nocturnal Fog

We have performed a one-dimensional and transient radiative heat transfer analysis in order to investigate interaction between atmospheric radiation and convective instability within a nocturnal fog. The radiation element method using the Ray Emission Model (REM²), which is a generalized numerical method, in conjunction with a line-by-line (LBL) method, is employed to attain high spectral resolution calculations for anisotropically scattering fog. The results show that the convective instability has a strong dependence on radiative properties of the fog. For the condition of a 20-m droplet diameter and liquid water content of 0.1 × 10^–3 kg m^–3;, the temperature profile within the fog becomes S shaped, and a convective instability layer forms in the middle or lower level of the fog. However, for the same water content and a 40-m diameter droplet, no strong convective instability layer forms, whereas for a 10-m diameter droplet a strong convective instability is observed.     

Perturbed physics ensemble using the MIROC5 coupled atmosphere–ocean GCM without flux corrections: experimental design and results

In this study, we constructed a perturbed physics ensemble (PPE) for the MIROC5 coupled atmosphere–ocean general circulation model (CGCM) to investigate the parametric uncertainty of climate sensitivity (CS). Previous studies of PPEs have mainly used the atmosphere-slab ocean models. A few PPE studies using a CGCM applied flux corrections, because perturbations in parameters can lead to large radiation imbalances at the top of the atmosphere and climate drifts. We developed a method to prevent climate drifts in PPE experiments using the MIROC5 CGCM without flux corrections. We simultaneously swept 10 parameters in atmosphere and surface schemes. The range of CS (estimated from our 35 ensemble members) was not wide (2.2–3.2?°C). The shortwave cloud feedback related to changes in middle-level cloud albedo dominated the variations in the total feedback. We found three performance metrics for the present climate simulations of middle-level cloud albedo, precipitation, and ENSO amplitude that systematically relate to the variations in shortwave cloud feedback in this PPE.     

Automated cloud tracking system for the Akatsuki Venus Climate Orbiter data

Japanese Venus Climate Orbiter, Akatsuki, is cruising to approach to Venus again although its first Venus orbital insertion (VOI) has been failed. At present, we focus on the next opportunity of VOI and the following scientific observations.We have constructed an automated cloud tracking system for processing data obtained by Akatsuki in the present study. In this system, correction of the pointing of the satellite is essentially important for improving accuracy of the cloud motion vectors derived using the cloud tracking. Attitude errors of the satellite are reduced by fitting an ellipse to limb of an imaged Venus disk. Next, longitude–latitude distributions of brightness (cloud patterns) are calculated to make it easy to derive the cloud motion vectors. The grid points are distributed at regular intervals in the longitude–latitude coordinate. After applying the solar zenith correction and a highpass filter to the derived longitude–latitude distributions of brightness, the cloud features are tracked using pairs of images. As a result, we obtain cloud motion vectors on longitude–latitude grid points equally spaced. These entire processes are pipelined and automated, and are applied to all data obtained by combinations of cameras and filters onboard Akatsuki. It is shown by several tests that the cloud motion vectors are determined with a sufficient accuracy. We expect that longitude–latitude data sets created by the automated cloud tracking system will contribute to the Venus meteorology.